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ABSTRACT 


This publication presents an antenna system that has been proposed 
as one of the candidates for the SCANSCAT (Scanned Scatter ometer ) 
radar application. It is the mechanically steered planar 
microstrip ref lectarray . Due to its thin, lightweight structure, 
the antenna's mechanical rotation will impose minimum angular 
momentum for the spacecraft. Since no power-dividing circuitry is 
needed for its many radiating microstrip patches, this electrically 
large array antenna demonstrates excellent power efficiency. In 
addition, this fairly new antenna concept can provide many 
significant advantages over a conventional parabolic reflector. 

The basic formulation for the radiation fields of the microstrip 
ref lectarray is presented. This formulation is based on the array 
theory augmented by the Uniform Geometrical Theory of Diffraction 
(UTD) . A computer code for analyzing the microstrip ref lectarray ' s 
performances, such as far-field patterns, efficiency, etc., is also 
listed in this report. It is proposed here that a breadboard unit 
of this microstrip reflectarray should be constructed and tested in 
the future to validate the calculated performance. The antenna 
concept presented here can also be applied in many other types of 
radars where a large array antenna is needed. 
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I . Executive Summary 


The mechanically steered flat microstrip reflectarray antenna 
has been proposed as one of the candidate antennas for the SCANSCAT 
(Scanned Scatterometer ) radar application. Due to its thin, light- 
weight structure, the antenna's mechanical rotation will impose 
minimum angular momentum for the spacecraft. Since no power 
dividing circuitry is needed for its many radiating microstrip 
patches, this electrically large array antenna demonstrates 
excellent power efficiency. As will be described in this report, 
this fairly new antenna concept can provide many significant 
advantages over a conventional parabolic reflector. 

The tasks that have been completed this year are the basic 
formulation of the antenna radiation fields and the generation of 
a computer code for analyzing the microstrip ref lectarray 's 
performances, such as far-field patterns, efficiency, etc. The 
formulation is based on the array theory augmented by the Uniform 
Geometrical Theory of Diffraction (UTD) . It is proposed here that 
a breadboard unit of this microstrip reflectarray should be 
constructed and tested in the future to validate the calculated 
performance. This antenna concept can also be applied in many 
other types of radars where a large array antenna is needed. 

II. Introduction 


1. Antenna Requirements: 

The Ku-band SCANSCAT antenna is required to generate two sets 
of pencil beams, as shown in Figure 1, to map the surface of the 
Earth for scatterometer application. One set should radiate at 36° 
from the nadir direction, while the other should be pointed at 49° 
from the nadir. Both sets of beams should be able to transmit and 
receive horizontally polarized electrical fields. Due to the 
azimuth scan speed of 16.4 rpm for the mechanical platform and 
relatively slower radar pulse return time, the transmit and receive 
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beams at 36° elevation need to be separated in azimuth by 0.59°. 
Similarly, the transmit and receive beams at 49° elevation should 
be separated by 0.77° in azimuth. Consequently, a total of four 
beams are to be generated at the frequency of 13.995 GHz. Each 
beam needs to have a 3-dB beamwidth of not more than 0.7° with a 
peak gain of 47 dBi or more and sidelobe level not higher than 
-20dB . 


2. Previously Proposed Concepts: 

A previously proposed antenna system was a mechanically steered 
dual-reflector configuration as illustrated in Figure 2. This 
configuration, although simple in design and relatively low in 
cost, will require large physical volume, mass, and, worst of all, 
a momentum compensation system that will result in even larger 
mass. As a consequence, a task was recently engaged to study other 
antenna concepts that will lead to smaller antenna volume, mass, 
and spin momentum. Quite a few concepts 111 have been investigated, 
which include electronically scanned phased arrays, hybrid 
electronically/mechanically scanned arrays, and mechanically 
steered arrays. The phased array approach was ruled out due to its 
extremely high cost and complexity. The only technically viable 
phased array approach is the active array with distributed T/R 
modules. With approximately 33,000 radiating elements, T/R 
modules, and phase shifters required, the whole antenna system will 
cost about $20 million. In addition to the high cost, it will be 
difficult (if not impossible) for the phased array to generate 
horizontally polarized beams in all the azimuth directions 
required. Also, an impractically large amount of D.C. power (in 
the order of many Kilowatts) may be needed to bias and control the 
large number of T/R modules and phase shifters. 

After a trade-off study, one antenna concept was selected for 
further investigation. This is the mechanically steered microstrip 
ref lectarray , which consists of a single flat array disc (3.1-meter 
diameter) that is space-fed by a small low-gain antenna and 
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mechanically spun in azimuth. This microstrip reflectarray offers 
much lower spin momentum, physical volume, and mass than the dual- 
reflector system. A detailed description of this antenna concept, 
as well as its advantages, are given in the following sections. 
Theoretical analysis of the antenna's performance was also carried 
out in this study effort. The basic array theory augmented by the 
Uniform Geometrical Theory of Dif fraction 12,31 was utilized for the 
theoretical analysis. A user oriented computer code for analyzing 
the microstrip reflectarray has been successfully used to produce 
many valuable data. 

III. Description of the Microstrip Reflectarray 

When many antenna elements with open or short circuited 
terminations are arranged in a planar aperture and are illuminated 
by a feed antenna as shown in Figure 3, these elements will re- 
radiate their illuminated energy into space. The total re-radiated 
energy will be non-coherent in phase, even when all the elements 
and their terminations are identical. This is because the fields 
that propagate to the elements from the feed have different path 
lengths, S lf S 2 , . ...S N , as shown in Figure 3(b), and thus form 
different phases. However, if each element's phase is adjusted to 
compensate for these different path lengths, the total re-radiated 
field can be made coherent and concentrated toward a specific 
direction. Many different types of radiators, such as horn, 
dipole, etc., can be used on the planar structure. If this planar 
structure is required to be thin, then a printed-circuit type of 
element needs to be employed. Several of these printed-circuit 
flat plate antennas have been attempted previously without much 
success. An early attempt 141 used the flat plate concept as a lens 
rather than as a reflector. Very low efficiency (25%) was 
reported. Recently, Malibu Research Center has demonstrated a flat 
plate reflector [51 without showing much data on efficiency. Both 
above attempts used either resonant slots or dipoles without any 
phase delaying transmission lines. Only the size of the slot or 
dipole is varied for phase trimming which is achieved by adding 
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reactance into the element's radiation impedance. This is 
equivalent to saying that the phase of the dipole's surface current 
is different for different dipole lengths. With properly designed 
phase distribution, beam coherence can be achieved by these arrays. 

This, however, will result in reduced radiation efficiency. Since, | 

for a particular frequency, there is only one optimum size of the 
resonant structure to transmit through or reflect energy, other 
sizes will result in low amplitude. Reference [4] printed that, j 

"amplitude must be sacrificed at the expense of phase change". A : 

more recent study [6! has used equal-size dipoles with periodic : 

spacing so that all the dipoles have the same illumination phase as 
the feed. This is the so-called frequency grating and is very i 

sensitive to frequency change. It is therefore used as a frequency | 

scanned offset-fed flat reflector antenna. Aperture efficiency of 
50% has been achieved by this antenna. Due to its required offset- 
fed configuration, this antenna cannot be effectively used in the 
SCANSCAT application. 

The flat plate antenna proposed here uses the ref lectarray 171 
concept where a short transmission line is attached to each 
radiator for phase adjustment. The electrical lengths of these 
transmission lines are made different depending on their radiators' 
positions from the feed. All the radiators, however, are identical 
in size and are made of thin microstrip patches. This antenna, as 
shown in Figure 4, is called the microstrip ref lectarray . It is 
composed of a thin (<0.02 wavelength) slab of dielectric material 
having one side completely covered with a layer of thin metal 
(which serves as a ground plane) and the other side etched with 
many identical metallic microstrip patches. A feed antenna, 
located at an optimally designed distance from the flat plate, will 
effectively illuminate all the metallic patches. The size of each 
patch, which can be rectangular, square, or circular, is made to 
resonate at the same frequency as the feed antenna. A short 
transmission line is connected to each patch at one end with the 
other end of the line either open or short circuited. This 
transmission line can be either a microstrip line etched on the 
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same side of the patches or a stripline sandwiched in an additional 
layered structure behind the ground plane. The advantage of the 
microstrip line is ease of fabrication with very little impact on 
antenna weight, while that of the stripline is minimum interference 
to the patch's radiation. When the radiation field of the feed 
antenna (in transmit mode) strikes each patch, the received 
resonant field of the patch will travel through its connected 
transmission line and be reflected by its open Or short circuited 
termination and then re-radiate through the patch into space. 
Thus, all the microstrip patches behave as re-radiators, while the 
short transmission lines serve as phase delay lines. The lengths 
of these transmission lines are intentionally made different for 
differently located patches so that the path delay differences from 
the feed antenna can be compensated. With proper design and 
calibration of these line lengths, the re-radiated fields from all 
the patches can be made coherent and concentrated toward the 
broadside direction. Also by re-designing the line lengths, the 
main beam can be directed toward other directions. Since the 
required phase changes for all the elements are between 0° and 
360°, the maximum length needed for the transmission line is only 
a half-wavelength. Consequently, the insertion loss associated 
with these short lines will be insignificantly small. However, 
this half-wavelength transmission line will work for a narrow 
bandwidth (<1%) application. For a wider bandwidth, a frequency 
excursion error will occur, especially for the outer elements of 
the array (assuming the feed is located at the center axis of the 
array) . In other words, the phase will accumulate more error for 
the outer elements. This accumulated phase error can be reduced by 
using longer transmission lines for the center elements and/or by 
using a larger f/D ratio (where f is the distance between the feed 
and the patches, and D is the diameter of the ref lectarray) . 

Since the microstrip reflectarray does not require any power 
divider, its efficiency in a large array system is much higher than 
a conventional array having the same aperture size. One possible 
drawback of this reflectarray is that, in addition to the re- 
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radiated fields from the patches, there will also be scattered 
field from the patches, reflected field from the ground plane 
(especially at off-resonant frequencies of the patch) , and 
diffracted fields from edges of the flat plate. These 
backscattered fields may increase the sidelobe level and possibly 
distort the main beam shape. However, as long as the aperture 
directivity of the flat plate is sufficiently higher (25dB or more) 
than the feed directivity, the backscattered energy will be 
insignificantly small. In other words, as will be demonstrated in 
a later section, the microstrip reflectarray will be an efficient 
antenna system only if it has a large number of elements (thousands 
or more) . 

IV. Application to SCANSCAT 

The proposed microstrip reflectarray for SCANSCAT is configured 
in Figure 5 where each patch has two orthogonally connected 
transmission lines with different lengths. Each one corresponds to 
a different elevation beam. Although not required, this orthogonal 
connection and the proper phase settings can separate the two 
different elevation beams 90° apart in azimuth. To meet the four- 
beam requirement of SCANSCAT, four feeds are needed as shown in 
Figure 5. Each two with the same polarization, due to their 
spatial separation, will generate the two separated transmit and 
receive beams. This is because each different feed location will 
form a different set of phases to illuminate the patches and, thus, 
generate a different beam position. Each two feeds that have 
different polarizations will generate two beams at different 
elevation and azimuth angles due to the two different settings of 
transmission line lengths on the microstrip patches. Because of 
the elegance of this design, all four beams will have horizontal 
polarization as required. These four beams are to be scanned in 
azimuth by spinning mechanically at the center of the circular flat 
plate. Due to the fact that the mass of the circular flat plate is 
more concentrated toward the center of rotation than the previously 
proposed dual-reflector design (Figure 2) , the spin momentum will 
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be much less. In addition, due to the microstrip design, the 
overall antenna volume and mass will also be much less than the 
dual-reflector system. 

To meet the SCANSCAT's 47dBi of minimum gain requirement, 
approximately 65,800 microstrip patches are needed. This huge 
number of elements does not impact the efficiency of the antenna 
system since no power-dividing circuitry is needed for the 
ref lectarray . The entire array has a diameter of 3.1 meters with 
a preliminarily designed feed location of 1.5 meters from the 
array. This 1.5 meters of feed separation (0.5 f/D ratio) is not 
necessarily an optimum number. If the feed gets too close to the 
array, the frequency excursion error will increase and therefore 
will narrow the operational bandwidth. If the feed is too far away 
from the array, the mechanical deployment and structure support 
problems may become severe. Certainly, the feed spillover 
efficiency and illumination efficiency are also determining factors 
in designing the feed location. This feed location optimization 
should be done together with a mechanical engineer. 

V. Analysis of the Microstrip Reflectarray 

Consider a planar array consisting of M x N microstrip patch 
elements that is non-uniformly illuminated by a low-gain feed at r f 
as shown in Figure 6. Let the desired beam direction be specified 
by unit value U Q . Then the re-radiated field in the U direction 
will be of the form 


M N 


E ( u) 'EE F(r mn * i f ) A(r n 


tf 0 > 


A(u 


S 0 ) 


/n*l n-1 


(i) 


exp {-jk[\r mn - r f ! + • u] + jaj + E z + E d 
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where F is the feed pattern function, A is the pattern function of 
the microstrip patch on the flat plate, r^ is the position vector 
of the ninth patch, and is the required transmission line phase 
delay of the mnth element for beam coherence. The condition that 
the beam will be coherent at desired direction U„ is 


a mn - k[\r mn - z f \ + r OT • u 0 ) = 2 nn, n = 0, 1, 2 


( 2 ) 


| 


The feed function F is modeled by cos q 0 function. For the pattern 
function A of the single square or rectangular microstrip patch on 
the flat plate, a simple closed form model using the dual-slot 
theory 131 is employed. This simple model, which is accurate enough 
for large array prediction, allows the computation time of many 
thousands of array elements to be significantly reduced. The | 

radiation patterns of the 0 and <p components from each slot of the 
dual-slot model, as illustrated in Figure 7, is given as follows: 

E sin (k a coscfr sin9) . cos ( k b sin<t> sin6) cos 0 

6 k a cos<|> sin0 (kb sin<|> sin0) 2 - (rc/2) 2 ( 3 ) 


E = - S3 - n US. a cos< t > sin0) 
* k a cos<{> sin 0 


cos (k b s in$_ sinQ) sin x cos0 

■ . > • ^ i *5 * 


(kb sin<J> sin0) 2 - (n/2) 


(4) 


where a is half the dielectric thickness (slot width) and b is half 
the patch width (slot length) , and K = 2n /77/X 0 . The terms E r and 
E d in equation (1) are, respectively, the specular reflected field 
from the flat ground plane and the diffracted field from the edges 
of the ground plane. Both E r and E d are calculated via the 
technique of the Uniform Geometrical Theory of Dif fraction 12,31 . 
Although E r does not give an accurate solution for the scattered 
fields from the patches and the ground plane, it gives the worst 
solution (maximum backscattered field is reflection) . This can be 
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assumed because the patches are separated within a very small 
distance (<0.02 wavelength) from the ground plane, and the flat 
plate reflectarray can thus be treated RF wise as a perfect 
conducting plane. This is especially true at off-resonant 
frequencies. To accurately predict the backscattered fields from 
thousands of non-uniformly illuminated microstrip patches with 
unequal lengths of microstrip transmission lines will require the 
development of a complex analysis technique which is beyond the 
scope of the current allocated fund. The fields E r and E d , no 
matter how accurate their calculation, are insignificant to the 
main beam when the array aperture gain is about 30 dB higher than 
the gain of the feed antenna, which is indeed the case for 
SCANSCAT. Nevertheless, these two terms are included in the 
analysis so that the worst possible sidelobe level due to scattered 
fields can be predicted. 

The efficiency of the microstrip reflectarray is primarily 
governed by the aperture illumination efficiency and feed spillover 
ef f iciency [8] . Aperture illumination efficiency is caused by unequal 
illumination of the array due to the feed pattern. The spillover 
efficiency is the ratio of the amount of feed energy that 
illuminates the entire array to the amount of energy that spills to 
the outside of the array. The calculation of these two 
efficiencies are very similar to that of a parabolic reflector. My 
colleague, Dr. Vahraz Jamnejad, with his many years of reflector 
experience, has assisted me in calculating these two efficiencies. 
With the q factor in the feed pattern cos q 0 chosen to be 3 and the 
feed separation of 1.5 meters from the 3.1-meter diameter array, 
the efficiency of the overall microstrip reflectarray for the 
SCANSCAT is calculated as following: 
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Table 1. Estimated microstrip reflectarray 
efficiency for SCANSCAT 


Type 

of Efficiency 

Efficiency in Percent 

spillover 

91 

illumination 

83 

patch loss 

95 

feed loss 

95 

termination loss 

95 

total 

65 


Because the proposed microstrip reflectarray for SCANSCAT is an 
electrically large array, the above estimated 65% antenna 
efficiency is considered quite good. In addition, the illumination 
and spillover efficiencies may be improved by optimally designing 
a special feed instead of using the cos q 9 feed pattern. 

■ 

A computer code based on the formulation in Equation (1) has 
been written with the FORTRAN language and user friendly inputs. I 

For the SCANSCAT antenna, the patterns of the 36° and 49° scanned § 

beams are calculated. For brevity's sake, only the 49° scanned i 

beam patterns are presented here. Figure 8 gives the coordinate : 

system of the reflectarray geometry. The radius of the | 

reflectarray is 72 X 0 and the patch elements are spaced 0.5 X 0 | 

apart, with X 0 being the free space wavelength. Each patch element § 

is 0.32 X 0 x 0.32 X 0 square with a substrate thickness of 0.01 X 0 and J 

a dielectric constant of 2.2. The feed has a symmetrical cos q 9 1 

, ... I 

power pattern with q = 3 and an illumination edge taper of -9dB. % 

Figure 9 is the calculated far field pattern in the X-Z plane when j 

all the phase delay transmission lines on the patches are set for * 

a beam scan of 49° in the x-z plane. The feed is located at (XF, | 

YF, ZF) = (0., 0., 0.) or (x, y, z) = (0., 0., 72 X c ) with an f/D | 

ratio of 0.5. Both the feed and the patches are polarized in the 
Y-direction so that, when the z-axis of the array is pointed at 
nadir and the beam is scanned in the x-z plane, a horizontally 
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polarized beam is achieved. This pattern shows a 3-dB beamwidth of 
0.75° and peak sidelobe level of -3 2dB. The beamwidth of 0.75° may 
be reduced down to the required 0.70° by a feed optimization 
program. It needs to illuminate the array more uniformly to 
increase the effective aperture and thus reduce the beamwidth. 
Figure 10 gives the 0-plane pattern (constant 0 = 49° plane) with 
a 3-dB beamwidth of 0.64° and a peak sidelobe of -31dB. The 
calculated directivity of this 49° scanned beam is 50.3dBi which, 
after subtracting 65% of efficiency loss, indicates an antenna 
peak gain of 48.4dBi. This is about 1.4dB above the required 
minimum gain. 

Figure 11 gives the 0-plane pattern when the feed is displaced 
0.4 X 0 in the +y-direction (XF, YF, ZF = 0., 0.4 X 0 , 0.), while all 
other parameters of the array are kept the same as those for 
Figures 9 and 10. The beam shows a shift of 0.34° in the -0 
direction. With a second identical feed located at YF = -0.4 X 0 , 
another beam with a 0.34° shift in the +0 direction can be formed. 
Separate transmit and receive beams can thus be achieved for the 
SCANSCAT. The shifted beam shown in Figure 11 has a coma lobe at 
-25dB level which is still acceptable to SCANSCAT. This coma lobe 
distortion is a consequence of phase errors that resulted from all 
the phase delay lines with settings designed for a focal feed 
rather than an off-focal feed. 

To demonstrate the effect of backscattered field components — 
E r and E d of Equation (1) — Figure 12 presents the pattern for the 
same antenna (diameter = 144 X Q ) as in Figure 9 except with the 
horizontal axis expanded and with less sampling accuracy in the 
pattern. This pattern does not include the effect of E r and E d and 
shows no far-out sidelobes above -60dB (mainbeam peak is normalized 
at OdB) . With this same antenna design. Figure 13 shows the 
pattern that does have the terms E r and E d included. Many sidelobes 
are near the -4 OdB level which is much lower than the SCANSCAT 
requirement. Now, let us look at a smaller reflectarray with 
diameter of 10 X Q , f/D ratio of 0.5, and feed pattern q factor of 
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3. The far field pattern of this antenna when scanned to 36° is 
shown in Figure 14 where E r and E d terms are not included. For the 
same antenna, the scan pattern with E r and E d included is given in 
Figure 15 where sidelobes of -17dB level are observed. Figures 12 
through 15 have demonstrated that the microstrip reflectarray 
antenna will not have serious high sidelobe problems caused by the 
backscattered fields when the aperture directivity is significantly 
higher (such as 30dB) than the feed directivity. In other words, 
the microstrip reflectarray is a more effective antenna when it is 
electrically large. Figure 16 plots the calculated antenna gain 
(estimated loss included) versus the array aperture diameter. It 
shows that, at the SCANSCAT frequency, the antenna gain drops much 
faster with reduction in aperture size when the diameter is less 
than one meter. 

VI. Advantages and Disadvantages of the Microstrip Reflectarray 

The numerous advantages of the microstrip reflectarray, in 
addition to its excellent application to SCANSCAT, are separately 
discussed below: 

1. The reflectarray, being in the form of a microstrip antenna, 
can be fabricated with a simple, low cost, and accurate 
etching process. Being flat, the microstrip reflectarray will 
be more cost effective than the parabolic reflector when 
manufactured. For example, the special molding process that 
is generally required for fabricating a paraboloid is not 
needed for a flat antenna. 

2. Due to the fact that no power divider is needed, the insertion 
loss of thousands of microstrip patches in the reflectarray 
will be the same as the insertion loss of a single patch and 
thus achieve good efficiency. The efficiency is much better 
than that of a conventional array with the same aperture size 
and is estimated to be comparable to that of a parabolic 
reflector (55-75% efficient) . 
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The main beam of the microstrip reflectarray can be fixed to 
point at large angles (up to 60°) from the broadside 
direction, while a parabolic reflector can only have limited 
scan (several beamwidths) . Phase shifters can be placed in 
the phase delay transmission lines for electronic beam 
scanning. 

4. Since the antenna is a flat structure, its mass is likely to 
be less than a curved parabolic reflector with equal aperture 
size. It can be more easily mounted onto the surface of a 
structure, such as a spacecraft's main body or a building, 
with less supporting structure volume and mass when compared 
to a parabolic reflector. As a possible application, this 
microstrip reflectarray can be made as the world's largest 
array antenna with relative ease in construction and low cost. 
This can be true because a flat reflectarray, where the RF 
loss is not a function of array size, can be constructed on 
flat land (aperture parallel to the land) with its feed 
mounted on top of a high tower. The size of the antenna is 
only limited by the tower height. A feed height of 1000 feet 
with an f/D ratio of 0.5 will result in an array with diameter 
of 2000 feet. This antenna size at 10 GHz can produce a 
pencil beam of 0.0035° in beamwidth and 95dBi of directivity. 

One major disadvantage of the microstrip reflectarray is the 
narrow bandwidth of the microstrip patch. It is certainly no match 
to the wide band property of the parabolic reflector (theoretically 
infinite bandwidth) . With conventional microstrip patch design, a 
maximum of 3% bandwidth can be achieved. With special design, such 
as a dual-stacked patch 191 , the bandwidth can be increased close to 
10%. The disadvantage of the narrow bandwidth of the microstrip 
reflectarray can be overcome somewhat, by utilizing multiple- 
frequency operation. This is made possible for the flat 
reflectarray by using a multiple-layer design as shown in Figure 17 
where the larger patches serve as ground planes for many of the 
smaller patches. A microstrip antenna having triple-frequency 
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capability with a three-layer design has been successfully 
demonstrated^ 1 . Other multiple-frequency designs are also 
possible, such as multiple-ring patches and interlaced different- 
size patches, shown in Figure 18. 

VII. Conclusion 

The flat plate microstrip reflectarray antenna has been 
analyzed by the array theory augmented by the Geometrical Theory of 
Diffraction. Antenna performances, such as radiation pattern, 
directivity, efficiency, etc., have demonstrated that a SCANSCAT 
mechanically steered antenna is feasible. Its radiation efficiency 
is comparable to a parabolic reflector antenna, while having many 
mechanical advantages over a parabolic reflector. Due to its low 
RF loss characteristic, this antenna can be used in many large beam 
scanning radar applications. 
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? 

NADIR 



EACH BEAMWIDTH s 0.70° 

T/R BEAM SEPARATION = 0.77° 


EACH BEAMWIDTH £ 0.70° 

T/R BEAM SEPARATION = 0.59° 


• PEAK GAIN = 47 dBi 

• POINTING KNOWLEDGE = 0.02° 

• POLARIZATION = HH 

• SIDELOBE LEVEL < -20 DB 


Figure 1. SCANSCAT antenna requirements. 





Figure 2. Mechanically steered dual-reflector configuration. 


mcJk wnmnm eust 


-17- 


PRECEDSNG PAGE BLANK NOT FILMED 


J/AA 




Figure 3. Planar array reflector antenna. 


FEED 

ANTENNA 



CIRCULAR 

MICROSTRIP 



Figure 4. Flat plate microstrip ref lectarray. 
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• TWO BEAMS IN DIFFERENT ELEVATIONS ARE SEPARATED 90° IN AZIMUTH WITH 
BOTH HAVING HH POLARIZATION 

• 65,800 MICROSTRIP PATCHES CAN BE FABRICATED WITH SIMPLE CIRCUIT 
ETCHING PROCESS 

• THE ARRAY REFLECTOR HAS A THICKNESS OF 0.032", THICKER HONEYCOMB - 
TYPE OF STRUCTURE IS NEEDED TO SUPPORT THE THIN FLAT ARRAY 

• REQUIRED ANTENNA SURFACE ACCURACY 5 0,030” 


Figure 5. Microstrip reflectarray for SCANSCAT. 



Figure 6. Coordinate system for array theory. 
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LOWER FREQUENCY 
PATCH 



Figure 17. Dual-frequency double-layer microstrip reflectarray 
concept . 


MULTIPLE-FREQUENCY MICROSTRIP PATCHES 



LOWER 

FREQUENCY 

RING 



(a) MULTIPLE-RING DESIGN 


(b) INTERLACED-ELEMENT DESIGN 


Figure 18. Multiple-frequency microstrip reflectarray designs 
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3313929*SCAN(1).SCAT(42) 


1 

C . . . 

2 

c . . . 

3 

c . . . 

4 

c . . . 

5 


t 


7 


6 


9 


10 


11 


12 


13 


14 


15 


16 

c 

17 

C 11 

18 


19 

99 

20 

C 

21 

Cl 2 

22 


23 

C 

24 

Cl 3 

25 


26 

C 

27 

Cl 4 

28 


29 

C 

30 

Cl 5 

31 


32 

C . . . 

33 

c . . . 

34 

c . . . 

35 

c... 

36 

c . . . 

37 

c 

38 

Cl 6 

39 


40 

C 

41 

C7 

42 


43 

C 

44 

Cl 7 

45 


46 

C 

47 

Cl 8 

48 


49 

C 

50 

Cl 9 

51 


52 

C 

53 

C26 

54 


55 


56 



THIS PROGRAM COMPUTES THE RADIATION PATTERN OF A LARGE 
MICRCSTRIP REFLECTARRAY WITH CIRCULAR APERTURE# ALL INPUT 
LENGTHS ARE IN WAVELENGTH. PROGRAMMED BY DR. JOHN HUANG 
AT JET PROPULSION LABORATORY# 12/12/1989. 

DIMENSION DBT(50C)#DBP ( 5 CO)# XAX <500 ) #XX <2 )# YY (2) 

COMPLEX CJ#CJJ#EXX#EYY#EZZ#FACP#EX#ET#EZ#ETH#EPH 
COMPLEX ETHR#£PHR#ETHD#EPHD 

COMMON/DD1/EPS#T#UT2#THIR#PHIR#THOR#PHOR#SSP 
C0MN0N/DD2/C J#PI #TPI#D PR 
C OMMON/ 003/ X PS# Y PS# FX#FY#FZ# IP 
C J=(G. #1 . ) 

PI=3. 14159265 
TPI=2.*PI 
DPR=1 8C./PI 
CJJ=C1.E-15#1.E-15) 

WRITE (6# 1 1 ) 

PQRMAT C * MICROSTRIP A g B DIMENS IONS=? * ) 

REAOC5#99) AX#BY 
F3RMATO 
WRITE<6#1 2) 

FORMAT C * SUBSTRATE DIELECTRIC CONSTANT & THICKNESS*?*) 

REAO(5#99)£PS#T 

WRITE (6#1 3 > 

FORMAT < * ELEMENT SPACING & ARRAY APERTURE RADIUS*?*) 

R E AD (5# 99) DX#DY#RA 
WRIT£(6#14) 

FORMAT ( * FEED LOCATION 8 FEED COSINE 0 FACTORS*?*) 

READ(5#99)FX#FY#FZ#QX#CY 

WRITE(6#15 ) 

FORMATt* FEED POLARIZATION# IF X IP=1# IF Y IP*2*) 
READC5#99)IP 

IF IP=1# BEAM NEEOS TO SCAN TO Y DIRECTION OR PHS*90 DEG FOR 
HH POLARIZATION# SCAN TO X DIRECTION OR PHS=0 OEG FOR VV 
P0LARIZATIQN####IP=2 BEAM NEEOS TO SCAN TO X DIRECTION OR PHS=0 
DEG FOR HH POLARIZATION# SCAN TO Y DIRECTION OR PHS*9G OEG 
FOR VV POLARIZATION. 

WRITE (6#1 6) 

FORMAT ( * REQUIRED BEAM SCAN ANGLE*?*) 

READ<5#99)THS#PHS 

WRITE<6#7) 

FORMAT ( * PLOT X-AXIS 8 Y-AXIS LENGTH IN INCH=?*) 

READ(5#99)XDIM#YDIM 

WRITE(6#1 7) 

FORMAT < * PATTERN CUT ANGLE# THETA IC=1### PHI IC = 2*> 

REAOC 5# 99) IC 
URITE(6#1 8) 

FORMAT < * LEFT BOUNO# RIGHT BOUND# INCREMENT OF PLOT*?* ) 

REAO(5#99)THL#THR#ITH#PHL#PHR#IPH 

WRITE<6#19) 

FORMAT ( * WRITE OUT REQUIREC ELEMENT PHASE# IF YES# IW=1*> 

REAO(5#99)IW 

WRITE (6#26) 

FORMAT ( * GTD GROUND PLANE EFFECT# IF YES# IGD*1*) 

READ(5#99)IGD 

THSR=THS/DPP 

phsr=phs/dpp 


2^ INTENTIONS? mm 
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57 


SAI=ATAN!RA/FZ) 




58 


C$Al=C0S(SAI) 




59 


CSA2=COS(SAI/2.) 




60 


EKl=gx*AL0G10CCSAI) 




61 


EK=EK1/2./AL0Gl0!CSA2) 




62 


E FFl = ! 1 .-CSA2**(2.*EK))/£K/TANCSAI/2.) 



63 


EFF=!2.*EK*1.)*EFF1*EFF1 




64 


EFS=1 CSA2**<2.*(2.*EKM.» 




65 


EFI=£FF/EFS 




66 

C 

WRITE(6*53)EFS*EFI*EFF 




67 

C58 

FORMAT!* SPILL-OVER* ILLUMINATION*- 

TOTAL 

5FF, = • * 3F9. 3) 


63 


Ql=QX-1 

- 


: 

69 


TQ1=2.*QX+1. 




70 


EFS-1 .“CSAI**TQ1 



r 

71 


E c 1=!1 .-CSAI**OX)/QX 



; 

72 


£F2=(1 .-CSAI**Q1)/Q1 



l 

73 


EF3=C1 .-CSAI**TQ1)/TQ1 



- 

74 


EFI=((£F1+EF2)**2.)/EF3 



=- 

75 


EFI=EFl/2./TAN(SAI)/TAN(SAI) 




76 


EFF=€FI *€F5 

* 


i 

77 


W RITE! 6*59) EFS*EFI*EFF 



| 

78 

59 

FORMAT ! * SPILL-OVER* ILLUMINATION* 

TOTAL 

EFF . =* * 3F9. 3) 

* 

79 


NX = I NT (RA*2. 0000001 /QX ) *1 

, , ' 


i 

80 


NY-INT!RA*2. 00000 01 / DY)*1 




SI 


NX1=NX-1 



z 

82 


NT1=NY-1 



i 

83 


HFX-FL0AT(NX-1 )/2.*0X 



| 

84 


HFY=FL0AT!NY-1)/2.*0Y 



i 

85 


IFdC.EO. 2)G0 TO 65 




86 


IA=ITH+1 



: 

87 


DAN=!THR-THL> /FLOAT! ITH) 



1 

88 


GO TO 66 




89 

>5 

IA=IPH+1 



1 

90 


OANa(PHR-PHL) /FLOAT! IPN) 



5 

■ 

91 

66 

DO 500 I J*1 * I A 



| 

92 


EXX=C JJ 



m 

93 


EYY=CJJ 



— 

94 


EZZ=C JJ 

* 


= 

95 


IF!IC.EQ.2)G0 TO 67 

- 


= 

96 


THPL=THL+!IJ-1)*0AN 



i 

97 


PMPL-PHL 



| 

93 


GO TO 63 




99 

67 

THPL=THL 



s 

100 


PHPL=PNL+(IJ-1)*0AN 



5 

101 

63 

ir=o 



B 

102 


1C OUNT =0 



1 

103 


THPR=ABS!ThPL)/OPR 



■ 

1 

104 


phpr=phpl/dpr 



I 

105 


IF(THPL.LT.O.)PHPR=PHPL/OPR+PI 



i 

E 

106 


CT=COS(THPR) 



■ 

107 


CP=CQS(PHPfl) 




108 


ST=SIN<THPR) 



- 

109 


SP=SlN(PHPR) 




110 

31 

YP=FL0AT!NY-IY-1 ) *OY-HFY 




111 


IX=Q 




112 


IY=IY+1 




113 

32 

XP=— FLOATINX— IX-1)*0X+HFX 
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114 

115 

116 
117 
113 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 
153 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 


IF(ABS<XP).LT.1.E-3)XP S 1.E“8 
IFCABSCYP) .LT.1 .£-3>YP=1 .E-8 
RP=SQRT(XP*XP+YP*YP> 

RAP=RA*0. 0000001 

IF(RP.GT.RAP)G0 TO 62 

I COUNT =IC0UNT+1 

SX=XP-FX 

SY=YP~FY 

SZ=-FZ 

S=SQRT(SX*SX+SY*SY+SZ*SZ> 

UX=SX/S 

UY=SY/S 

uz=sz/s 

SXY=SQRTCSX*SX+SY*SY) 

THF=ACQS(FZ/S) 

PHF=ATAN2 (UYr UX) 

IF<PHF.LT.O.)PHF=PHF+TPI 

EPAT=CaS<PHF)*COS<PHF)*<COSCTHF))**QX+SIN<PHF)*SIN<PHF>* 

&<COSCTHF))**QY 

SS=SQRT<SZ*SZ+RP*RP) 

SFAC=SS-FZ 

SFA=FLOATCINTC$FAC)) 

PHAS1 =TPI*C5FAC— SFA) 

PHN=ATAN2<YP,XP) 

IF(PHN.LT.O.)PHN=PHN+TPI 
PHAS2=TPI*RP*COS(PH$R-PHN) *SIN(THSR) 

PHASR=PHAS2-PHAS1 

TES=PHASR/TPI 

ITES=INT(TES> 

PMASO=(TES“FLOAT CITES) )*TPI*OPR 
IFCIJ.EQ.1 .AN0.IW.EQ.1 >WRITE(6,*)XP,YP*PHAS0 
C ABOVE ARE PATCH LOCATION ANO REQUIRED PHASE OELAY 
IF ( IP. EQ. 2)G0 TO 61 
HT2=BY+T 
00 36 1=1 rZ 
XPS=XP-(AX+T>/2. 

IF(I„EQ.2)XPS=XPMAX + T)/2. 

YPS=YP 

XS=FX-XPS 

YS=FY-YPS 

ZS=FZ 

SSP=SQRT(XS*XS+YS*YS+ZS*Z$) 

UXS=XS/SSP 

UYS=YS/SSP 

UZS=ZS/SSP 

PHIR=ATAN2(-UYS,-UXS)+PI 
IFCPHIR.LT.O.) PHIR=PHI r+tp I 
THIR=AC0SCUZS) 

RSP=SQRT(XS*XS+YS*YS) 

PHNN=ATAN2(-UYS,-UX$) 

IF(PHNN.LT.G.)PHNN=PHNN+TPI 

PHASO=TPI*RSP*CQS(PHPR-PHNN)*SIN<THPR) 

FACP=CEXP(CJ*(PHASO-PHASR))*EPAT 

THOR=THPR 

PHOR=PHPR 

CALL SOURCEC£XP,EYP,EZP) 

CALL SLOTX(EXP/-EYP#EZP#EX,EY*EZ) 
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171 


EXX=EXX+£X*FACP 

172 


EYY=EYY+£Y*FACP 

173 


EZZ=EZZ*EZ*FACP 

174 

36 

CONTINUE 

175 


GO TO 62 

176 

61 

WTZ=AX+T 

177 


00 37 1=1,2 

178 


XPS=XP 

179 


YPS=YP-<BY+T)/2. 

18Q 


IFCI.EQ.2) YPS=YP+(0Y+T>/2. 

181 


XS=FX-XPS 

182 


YS=FY-YPS 

183 


IS-Fl 

184 


SSP=SQRT(XS*XS+YS*YS+ZS*ZS) 

185 


UXS=XS/$SP 

186 


UY$=YS/$SP 

187 


UZS=ZS/SSP 

188 


PHIR=ATAN2(-UYSx-UXS)+PI/2. 

189 


IF (PHI R .LT . 0 . ) PH IR=P HIR+TPI 

190 


THIR= ACOS(UZS) 

191 


RSP=$QRT<XS*XS+YS*Y$) 

192 


PHNN=ATAN2(-UYSx-UXS) 

193 


IF(PHNN.LT.O.)PHNN=PHNN+TPI 

194 


PHAS0=TPI*R5P*C0S< PHPR— PHNN) *SINCTHPR) 

195 


FACP=CEXP(CJ*<PHASO— PHASR) )*£PAT 

196 


TH0R=THPR 

197 


phor=pmpr-pi/2. 

198 


if(phor.lt.o.)phor=phor+tpi 

199 


CALL SOURCECEXPxEYPxEZP) 

200 


CALL SLOTX(EXPxEYPxEZPxEXxEYxEZ) 

201 


EXX=EXX-EY*FACP 

202 


EYY=E YY+EX*F ACP 

203 


EZZ=EZZ+EZ*FACP 

204 

37 

CONTINUE 

205 

62 

IX=IX+1 

206 


IF(IY.GT.NY)GO TO 31 

207 


IF<IX.GT.NX1)G0 TO 31 

208 


GO TO 32 

209 

81 

ETH=EXX*CP*CT*EYY*SP*CT-EZZ*ST 

210 


EPH=-EXX*SP+EYY*CP 

211 


IFIIGD.GT.1 )G0 TO 82 

212 


IF<IC.E0.2)G0 TO 82 

213 


IF<ABS<PHS) .GT.0.01 )G0 TO 82 

214 


IF(A6S<THPL) .LT.2.)G0 TO .32 

215 


ETHR=CJJ 

216 


EPHR=CJJ 

217 


RO=«=Z*TAN(THPR) 

218 


IF(RD.GT.RA)GO TO 83 

219 


XPS=RD 

220 


IF(THPL.LT.O.)XPS=-RO 

221 


YPS=0. 

222 


SSP=FZ/COS(THPR) 

223 


FAC=C£XP<-CJ*TPI*SSP) 

224 


EPAT=COS(THPR)**QX 

225 


PHASO=TPI*XPS*SIN(THPR> 

226 


IF(THPL.LT.C.)PHASQ=-PHASO 

227 


FACP=CEXP(CJ*PHASO)*EPAT*f AC 
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228 

229 

230 

231 

232 

233 

234 

235 

236 

237 
233 

239 

240 

241 

242 

243 

244 
24 5 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 
283 


CALL SOURCECEXP,EYP,EZP) 

EXX=-EXP 
E YY=“EYP 
EZZ=EZP 

ETHR=C£XX*CP*CT+EYY*5P*CT-EZZ*ST)*F ACP 
EPHR=<-5XX*SP+EYY*CP)*FACP 
83 TH0R=THPR 
PH0R=PHPR 

CALL E0GE(RA,QX,ETHD,EPH0) 

ETH=ETH+ETHR+ETHO 
EPM=EPH+EPHR+EPHO 
32 AET*C A8SCETH) 

AEP=CA8SCEPH> 

IFCAET.LT.1.£-15)AET=1 .£-15 
IFCAEP.LT.1 ,£-15>AEP*1 .E-1 5 
DBTCIJ)*20.*AL0G10(A£T) 

OaPCIJ>*20.*ALOGlGCAEP) 

XAX(IJ)=THPL 

IF(IC.EQ.2)XAXCIJ)*PHPL 

WRITE(6,77>XAXCIJ>,D8T(IJ),08PCIJ) 

77 FORMAT ( • ANGLE**,F12.5,3X, *DBT=*,F12.5,3X, , 0BP**,F12.5) 

500 CONTINUE 

WRITEC6,97)IC0UNT 

97 FORMAT!* TOTAL NUMBER CF ELEMENTS* *,I5) 

YMT=-1000. 

YMP=-1 000. 

00 22 1=1, IA 

IFCYMT.LT.08TCI))YMT*0BT (I) 

IFCYMP.LT.08PC I) )YMP=OBP<I) 

22 CONTINUE 
YMX=YMT 

IFCYMT.LT. YMP) YMX=YHP 

DO 21 1*1 , I A 

08TCI)=0BTCI)-YMX 

08PCI)*0BPCI)-YMX 

IFCOBTCI) .LT.-60.)OBTCI)=-60. 

IFC0BPCI).LT.-60.)0BPCI)*-60. 

21 CONTINUE 

CALL BGNPLT 

CALL PLFQRMC*LINLIN*,X0IM,YDIM) 

XX Cl ) *THL 
XX(2)*THR 

IFCIC.EQ.2)XXC1)=PHL 

IFCIC.£Q.2)XXC2>*PMR 

YYC1)*-60. 

YY C 2) =0. 

CALL PLSC ALCXX, 2, 040504, YY *2,060506) 

CALL PLGRAF(*MICROSTRIP REFLECTARRAY *,* ANGLE* , *08* ) 

CALL PLCURVCXAX,0BP,IA,0,0) 

CALL PLNUP 
CALL PLNTYPC5) 

CALL PLNDNCO.,0.) 

CALL PLCURVCXAX,08T,IA,0,0) 

CALL ENOPLT 

STOP 

ENO 
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331 3929*SCAN<1). SLOT X(5) 

1 SUBROUTINE SLOTX(EXP*EYP*E ZP*EX#EY* EZ) 

2 C...THIS SUBROUTINE GIVES 3£”R AOIATED FIELD FROM ONE EOGE SLOT 

3 C...OF A MICROSTRIP PATCH 


4 


COMPLEX Cj,FAC,EXrEY,EI 

5 


COHMON/OD1/EPS<T,HL,THIR,PHIR,THOR,PHOR,RI 

6 


CQHH0N/00 2/CJ,t>I,TPI#0PR 

7 


SP$=SQRT(EPS> 

3 


A=T/2. 

9 


3-HL/2 • 

10 


CPHI=COS(PHIR) 

11 


SPHI=SINCPHIR> 

12 


CTHI-CQS(THIR) 

13 


STHI=SIN<THIR> 

14 


ARG1=TPI*A*CPHI*STHI*SPS 

15 


ARG2=TPI*B*SPHI*STHI*SP$ 

16 


IFURG1 .LT.1 .E-4JG0 TO 11 

17 


F1=SIN<ARG1)/ARG1 

18 


GO TO 12 

19 

11 

FI *1 • - 1 

20 

12 

IF ( ABS( ARG2-PI/ 2. ) .LT. 1 . E-4)G0 TO 13 

21 


F2=C0S<A«G2)/CARG2*ARG2-PI*PI/4.) 

22 


GO TO 14 

23 

13 

F2=-1 ./PI 

24 

14 

EPH=F1*F2*SPHI*CTHI 

25 


£TH=“F1 *F2*CPHI 

26 


EXT=ETH*CTHI*CPHI-£PH*SPHI 

27 


EYT=ETH*CTHI*SPHI*EPH*CPHI 

23 


E2T=-ETH*STHI 

29 

C 

EII=EXP*EXT 

30 


£II=EXP*EXT+EYP*EYT+EZP*EZT 

31 


CPH=COS(PH0R) 

32 


SPH=SIN(PHOR> 

33 


CTH=C0S(TH0R) 

34 


STH=SINCTH0R) 

35 


ARG1=TPI*A*CPH*STH*SPS 

36 


arg2=tpi*b*sph*sth*sps 

37 


IFIARG1 .LT.1 .E-4)G0 TO 17 

38 


F1=SIN<A8G1)/ARG1 

39 


GO TO 13 

40 

17 

F 1 =1 . 

41 

18 

IF { ABS< ARG2-PI / 2. ) .LT. 1 .E-4) GO TO 19 

42 


F2=C0S< ARG2)/(ARG2*ARG2-PI*PI/4.) 

43 


GO TO 20 

44 

19 

F2=-1./PI 

45 

20 

EPH=EII*F1 *F2*SPH*CTH 

46 


£TH=£II*(-F1 *F2*CPH) 

47 


FAC=CEXP< -CJ*TPi*RI)/ R I 

48 


EX=(ETH*CTH*CPH-EPH*SPH)*FAC 

49 


EY=(ETH*CTH*SPH+EPH*CPH)*FAC 

50 


Ez=-eth*sth*fac 

51 


RETURN 

52 


END 


SPRTsS TT. SOURCE 
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3313929*SCAN<1).S0URCE<1> 

1 SUBROUTINE SQURCE{£XP*EYP*EZP) 

2 C...THIS SUBROUTINE GIVES SOURCE FIELD IN PATCH COORDINATE SYSTEM 


3 


C0MNQN/003/XP* YPrFX^FY rFl* 

A 


IF (IP. EQ. 1 ) GO TO 11 

5 


SX=XP-FX 

6 


SY=FZ 

7 


SZ=YP— FY 

S 


GO TO 12 

9 

11 

S X=FZ 

10 


SY=YP“FY 

11 


SZ=XP-FX 

12 

12 

SS=SQRT(SX*SX+SY*SY*SZ*SZ) 

13 


UXS=SX/SS 

14 


UYS=SY/SS 

15 


UZS=SZ/SS 

16 


THR=ACOS(UZS) 

17 


PHR=ATAN2(UYS*UXS) 

18 


IF<PHR.LT.O.)PHR=PHR+TPI 

19 


EX=— COS (THR) *COSCPHR > 

20 


EY=-COS<THR)*SIN<PHR) 

21 


EZ=SIN(THR) 

22 


IFCIP.EQ.1 )G0 TO 14 

23 


EXP=EZ 

24 


EYP*-EX 

25 


£ZP*-EY 

26 


GO TO 15 

27 

14 

£XP=EZ 

28 


EYP=EY 

29 


EZP=-EX 

30 

15 

RETURN 

31 


ENO 

3PRT,S 

TT.EDGE 
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12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


EDGE (4) 

SUBROUTINE EOGE(RA,QX, EPNO^EPRO) 

SINGLE EOGE DIFFRACTED FIELDS 
COMPLEX CJ,EPNO*EPRO,FACP,FAC 
COMPLEX OSxOH^DPSxOPH 

COMMON/DD1/£PS,T,HL/-THIR,PHIR,THOR,PHOR,RI 

COMMON/002/CJrPI/TPIxOPR 

COMMON/ OOZ/XP,'IP,fX,FX ,FZ, IP 

XP=RA 

YP=0. 

CALL SOURCECEX^EY^EZ) 

THOR=ATAN(FZ/RA) 

EPAT=C0SCPI/2.-TH0R)**QX 

EPNI=EX*SIN(THDR)*EZ*COS(THOR) 

EPRI-EY 

R=RA/COSCTHDR> 

FAC=CEXP(-CJ*TPI*R)/R*EP AT 
PHP=THDR*OPR 
PH=TH0R*0PR+90, ' '' 

IFIPHOR.GT.O.OI )PH=9Q.-TH0R*QPR 

CALL OU 10S,DH, OPS,OPH,a, PH,PH? ,90 . ,2 . ) 

ENI=COSCTHOR) 

ENS=C0S(Pl/2.-TH0R) 

IFCPHOR.GT.O.OI )cNS=C0S(PI/2.+TH0R) 

SGD=C1 ,/R)-(CENI-ENS>/RA) 

RH0-A8S(1./SG0) 

FACP=C£XP<CJ*TPI*RA*SIN(TMCR>) 

IF<PHOR.GT.0.01)FACP=CEXP(-CJ*TPI*RA*SIN<THOR)) 

EPND=EPNI*PAC*QH*SQRT ( RHO) *F ACP 

EPRD=EPRI*FAC*OS*SQRT<RHO>*FACP 

IF(SGD.LT.O.)EPND=EPND*CJ 

IF(S&O.LT.O.)£PRD=£PRO*CJ 

XP=-RA 

CALL $OURC£(EX,EY,£Z> 
EPNI=EX*SIN(THOR>-£Z*COS(THOR) 

EPRI=EY 

PH*90.-TH0R*0PR 

IF<PH0R.GT .0. 01) PH=THOR*DPR+90. 

CALL DW<OS*OH,OP$,OPH*R,PH*PHP*90.*2. ) 
ENS=C0S(PI/2.+TH0R) 

IF { PHOR.GT .Q.01)ENS=CQ$(Pl/2.-THQR) 

SG0=<1 ,/R)-<(ENI-ENS)/RA> 

RH0=A3S(1 ,/SGD) 

F A CP=C EXP (“CJ*RA*TPI *S IN (THOR) ) 

Ir( PHOR.GT. 0.01 > F ACP=C £X P( CJ*RA*TPI *S TN (T HOR ) ) 

EPN0=EPN0+6PNI*FAC*0H*S0RT(RHO)i*FACP 

E PRO=EPRD+£PRI*FAOOS* SORT <RHO>*F ACP 

IF(SGD.LT.O.)EPNO=EPNO*CJ 

IF(SGO.LT.O.)EPRO=EPRO*CJ 

RETURN 

END 
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3313929 *SCAN(1).DW<0) 

1 SUBROUTIN' OW ( DSr DHs DPS#DPHrR*PH/’PHP*BO*FN) 

2 C *** WEDGE DIFFRACTION AND SLOPE DIFFRACTION COEFFICIENT *** 

3 C *** FOR THE SOFT AND HARD E.C. *** 

4 COMPLEX D I N^ DI P^ DPN# DP P * DS*DH*DP5*D PH 

5 BETN=PH-PHP 

6 CALL DI(DIN*R*BETN*BDrFN) 

7 CALL DPI(DPN*R^BETN^BD*FN) 

8 I*=(ABS(PHP) .GT.2.5E-4. AND. ABS(PHP).LT. (FN*1 80.-2. 5E-4) 

9 8) GO TO 10 


10 


DS=(0.#Q.) 

11 


DH=OIN 

12 


DP$=DPN 

13 


DPH=<0.#0.) 

14 


RETURN 

15 

10 

CONTINUE 

16 


BtTP=PH+PHP 

17 


CALL DICDIP/R*BETP*80#FN) 

18 


CALL 0PI<DPP#R*6£TP/’BQrFN> 

19 


DS=DIN-DIP 

2C 


DH=DI N+DIP 

21 


DPS=OPN+DPP 

22 


OPH=OPN-OPP 

23 


RETURN 

24 


ENO 
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22 

23 

24 
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26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 
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43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 


C *** 

c 


442 

443 


542 

123 


I<2> 

SUBROUTINE 01 (DIR,R,-BET,.BO,FN) 

INCIDENT (BET=PH-PHP) OR REFLECTED (BET=PH+PHP) *** 
PART OF WEDGE DIFFRACTION COEFFICIENT *** 

COMPLEX TOP*COM,EX*UPP I^UNPI * FA* DIR 

DATA PI,TPI,DPR/3.14159265,6.2831853*57.29577958/ 

ANG-BET/OPR 

T0P=-CEXP(CMPLX(0.*-PI/4.)> 

0EM=2.*TPI+FN*SIN(80/DPR) 

COM=TOP/DEH 

SQR*SQRT(TPI*R> 

DNS=(PI+ANG)/(2.0*FN*PI) 

SGN=SIGN(1./.DNS> 

N=IFIX CABS CONS) +0.5) 

DN=SGN*N 

A= 1 . 0+C 0 S ( A NG- 2 . 0 *FN * P I * ON ) 

BOTL * 2.0*SQRT(ABS(R*A) ) 
EX=CEXP{CHPLXCO,0*TPI*R*A>) 

call frnels itVSVKoTU 

C=SQRT(PI/2.0)*{0.5-C) 

S- SORT{PI/2.0)*(S-0.5) 

FA=CmpLX(0.,2.)*SQR*6X*CMPLX(C*-S> 

RAG=<PI+ANG)/(2.0*FN) 

TSIN=SIN(RAG) 

TS=ABSCTSIN) 

IFCTS.GT.1 .E-5) GO TO 442 
COTA=-SQRT(2.0)*FN*SlN(ANG/2.0-FN*PI+DN) 

IF (COSC ANG/2 •0“FN*PI+0N) .LT.O.C) COTA=-COTA 
GO TO 443 

COTA=SQRT<A)*CO$(RAG)/TSIN 

UPPI=COM*COTA*FA 

0NS=(-PI+ANG)/(2.0*FN*PI) 

SGN*SiGN(1.,DNS> 

N=IFIX(ABS CONS) +0.5) 

ON-SGN*N 

A=1 . Q+COSC ANG~2 . 0*FN*PI*DN) 

BOTL * 2.0+SQRTCABSCR+A)) 
EX=CEXP(CMPLX<0.0,TPI*R*A>) 

CALL FRNELS CC*$*BOTL> 

C = SORT (PI/2.0) *(0.5-0 
S= SORT (PI/2 • 0) *( S“0. 5 ) 
FA=CMPLX(0.,2.)*SQR*EX*CHPLX(OS) 
RAG=(PI“ANG)/C2.0*FN) 

TSIN=SIN(RAG) 

T S-ABS(TSIN) 

IF(TS.GT.1.E-5) GO TO 542 
COT A= SQRT(2.0)*Fn*SIN(ANG/2.0-FN*PI*DN) 
IF(COS(ANG/2.0-FN+PI*ON) .LT.0.0) COTA=-COTA 
GO TO 123 

COTA=SORT(A)*COS( RAG)/TSIN 

UNPI=COM*CQTA*FA 

OIR=UPPI+UNPI 

RETURN 

END 
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3313929 

*SCAN<1).0PI (0) 

1 


SUBROUTINE DPI<DPIR*R*BET*80*FN) 

2 

c *** 

INCIDENT <BET=PM-PHP) OR REFLECTED <BET=PH+PHP) *** 

3 

c *** 

PART OF riEDGE SLOPE DIFFRACTION COEFFICIENT *** 

4 


COMPLEX TOP*COH*EX*UPPI*UNPI*FPA*OPIR 

5 


OATA PI*TPI*DPR/3.14159265*6.2831853*57.29577958/ 

6 


ANG*8ET/0PR 

7 


SB0=SIN(80/0PR> 

8 


T0P=CEXP<CMPLX(0.*-PI/4. )) 

9 


0EN=4.*TPI*f : N*FN*SB0*SB0 

10 


C0M=T0P/0EM 

11 


DNS- < PI+AN6) / ( 2 • 0*FN*P I) 

12 


SGN=S IGN< 1 . *0NS) 

13 


N-IFIX(ABS£DNS)*0.5) 

1* 


ON s SGN*N 

15 


A=1.0+C0S<ANG-2.0*FN*PI*GN) 

16 


BOTL * 2.0*SQRTiABSfR*Ar) 

17 


EX=CrXP<CMPLX<O.Q*TPI*R*A) ) 

18 


CALL FRNELS CC*S*90TL) 

19 


C = SQRT< PI/ 2.0) *<0.5-0 

20 


S= SORT (PI/2.0)*(S“0.5) 

21 


FPA=TPI*R*(CMPLX(0.*2. )+4. *SQRT<ABS<TPI*R*A))*EX*CMPLX<C*S)) 

22 


RAG=<PI+ANG)/C2.0*FN) 

23 


TSIN=SIN(RAG) 

24 


T S = TSIN*T SIN 

25 


IFCTS.GT.1 .E-5) GO TO 442 

26 


C SC A=— 2 . *FN*FN*COS< ANG-TPI *FN*DN) /COS ( ( PI +ANG) /FN) 

27 


GO TO 443 

28 

442 

CSCA=A/TS 

29 

443 

UPPI=COM*CSCA*FPA 

30 


ONS=<-PI+ANG)/<2.0*FN*PI) 

31 


SGN=S IGNC1 • *DNS) 

32 


N=IFIX<A8S<DNS)+0.5) 

33 


ON s SGH*N 

34 


A=1.0+COS(ANG-2.0*FN*PI*CN) 

35 


BOTL s 2.0*SQRTCABS(R*A)) 

36 


£X=CEXP<CHPLX<Q.O*TPI*R*A)) 

37 


CALL FRNELS <C*S*BOTL) 

38 


C = SQRT<PI/2.0) *<0.5-0 

39 


S* SQRT <PI/2.Q)*<S-Q.5 ) 

40 


PPA=TPI*R*(CHPLX<0.*2.)+4. *SQRT<ABS<TPI *R *A) ) *EX*CMPLX <C*S ) ) 

41 


R AG=( PI-ANG) / <2 .0*FN ) 

42 


TSIN*SIN(RAG) 

43 


TS=TSIN*TSIN 

44 


IF<TS.GT.1 .E— 5> GO TO 542 

45 


CSCA=-2.*FN*FN*C0S<ANG-TPI*FN*0N) /COS << PI-ANG) /FN) 

46 


GO TO 123 

47 

542 

CSCA=A/TS 

48 

123 

UNPI=COM*CSCA*FPA 

49 


OPIR=UPPI-UNPI 

50 


RETURN 

51 


ENO 

SPRT*S 

TT.FRNELS 
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51 3929*SCAN<1 ) .FRNELS<0) 

1 SUBROUTINE FRNELS (C # S# X S ) 

2 C THIS IS THE FRESNEL INTEGRAL SUBROUTINE WHERE THE INTEGRAL IS FROM 

3 C U=Q TO XS# THE INTEGRANO IS EX?<-J*PI/2.*U*U)#ANO THE OUTPUT IS 

4 C C<XS)-J*SCXS) . 

5 DIMENSION A ( 1 2 ) # B ( 1 2) # CC <1 2) ,0 C 1 2 ) 

6 OAT A A / 1. 595769 1 40#-0. C00001 702#-6. 80 35 68 554,-0. 000 5763 61 * 6 . 9 20691 

7 *902# -0.01 68986 57 #-3.050485660# -0.07 575241 9# 0 . 850663781 #-0.02563904 

8 *1#-O.150230960#0. 034404779/ 

9 DATA B/-0. 000000033#4. 25 5387 524#-0. 3000 92 31 0#-7. 7800 2C400# -0.0095 2 

10 *0395# 5.0751 6 1298#-0.1 38341 947#-1. 3637291 2 4#-0. 40334? 276#0. 70222201 

11 *6#-0.216195929#0. 019547031/ 

12 DATA CC /Q. #-0. 02493397 5 #C. C00003936#0. 005 7709S6#0. 00068989 2#-0. 009 

13 *4971 36#0. 01 1948809#-0.306748873#0. 0002464 20#0. 0021 02967#-0. 001 2179 

14 *30#0. 000233939/ 

1 5 DATA 0/0. 1 99471 140#0.0 000000 22#-0. 0093 5 1341 #0. 0000 23006 #0.00485 146 

16 * 6 # 0.001 90321 8#-0. 0171 2 29 1 4#0 -029064067#-0 .02792395 5# 0 . 01 649 7303 #-0 • 

17 *.005598515#0. 000838386/ 

18 DATA PI /3. 1 41 59265/ 

19 IFCXS.LE.O.O) GO TO 414 

20 X = XS 

21 X * PI*X*X/2.0 

22 FR=0.0 

23 FI=0.0 


24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
33 

39 

40 

41 

42 

43 

44 

45 

46 

47 
43 
49 


A = 13 

IF(X-4.0) 1 Q#4Q#4Q 
10 Y=X/4.0 

20 K=K-1 

FR=<FR+A(K))*Y 
FI=CFI+BOO)*Y 
IFCK-2) 30#30#20 
30 FR=FR+A<1) 

FI=FI*B<1 ) 

C=(FR*COS(X)+FI*SIN(X) >*SQRT(Y) 
S=<FR*SIN<X)-FI*C3S(X))*SQRTCY) 
RETURN 

40 Y=4.0/X 

50 K=K-1 

FR=(FR+CC(K))*Y 
FI=<FI+D<K>)*Y 
IF(X-2) o0#60#53 
60 FR=FR+CC(1) 

F I =FI+0 (1 ) 

C=Q.5+<FR*C0S(X)+FI*SIN(X) )*SQRT(Y) 
S=0. 5+<FR*SINCX)-FI*C0S(X) >*SQRt(Y) 
RETURN 
414 C=-0.0 

S=-0.0 
RETURN 
END 


TP=1 0. 363 SUP=6 . 1 07 CPU=.009 10=2.609 CC-ER=3.488 
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